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Cox-Thompson fixed-energy quantum inverse scattering method is developed further to treat long 
range Coulomb interaction. Depending on the reference potentials chosen, two methods have been 
formulated which produce inverse potentials with singular or finite value at the origin. Based on the 
quality of reproduction of input experimental phase shifts, it is guessed that the p — a interaction 
possesses an interesting repulsive hard core. 
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1 Introduction 

Inverse quantum scattering methods [TJ [3J |31 H] represent useful techniques to assess effective interactions 
between colliding composite particles in a model independent way. In nuclear physics the modified 
Newton-Sabatier (mNS) method [SJ [3] has been used to determine, in general, complex valued optical 
potentials describing various systems, e.g., "'^^C - ^^C elastic scattering [71 H] or to guess the spin-orbit 
potentials arising in the p - a and n - a collisions [5]. The Newton-Sabatier (NS) method [IHl [EL has 
the property that it requires an infinite set of phase shifts otherwise the first moment of the potential 
generated vanishes. This induces an unpleasant oscillation of the potential at large distances, in addition 
to the artificial singularity at the origin which remains to be present in the potentials generated also by 
the mNS method. 

In the recent decade another procedure, the Cox-Thompson (CT) inverse quantum scattering method 
[12] has been investigated [ISl dH US] • This method has the advantage that it requires a finite set of phase 
shifts and possesses a non-zero first momentum of the potential generated. In addition to this, in general, 
it also reproduces test potentials better than the mNS method [2] ■ The CT method has been applied so 
far to uncharged particle scattering in order to reconstruct short ranged potentials from synthetic phase 
shifts 13 , and also to construct n - a potential from experimental data [14j . We shall develop here 
the method further with the intention that the CT method can be applied to charged particle scattering 
when the potential to be constructed has a long ranged Coulomb tail. The various extensions will be 
applied to both synthetic and experimental phase shifts. The results are obtained by using the nonlinear 
solver package of MATLAB. 



The original CT method belongs to the fixed-energy procedures when a finite set of N phase shifts 
{Si}ieS is converted to the spherical potential V{r) appearing in the ^— th partial wave radial Schrodinger 
equation 
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2 Cox-Thompson (CT) method 




(1) 



with boundary conditions \I'z(fc;0) oc r'"*"^ and ^i{k;r — >■ oo) oc sin(fcr — + Si). 
We shall also refer to another form of the Schrodinger equation 

+ 1 - q{x)^ M^) - 1(1 + mi^l (2) 

where dimensionless quantities have been introduced for the distance x — kr, for the potential q{x) — 
E^^V{xk^^) with E — ?i^fc^/(2m), and the dimensionless wave function ipi{x) is related to the wave 
function ^';(fc;r) through the asymptotic normalization constant \E'i(fc;r) = Ci^i{kr). 

The conversion from phase shifts {(5;} to potential V{r) can be carried out in two steps [131 [T3]. First, 
we solve either of the following two systems of highly non-linear equations 

i-./cr 2 + z(/c+-/cr)' ^' 

with 

/C±= Yl [Msin]^L[A/-i]Lpe±'('^'>/^ I e S. (4) 

and 



Mcos L{L + 1) - + 1) \ cos {{I - i)f ) 

for the set T of the shifted angular momenta {L} where the relations 5 n T = and |r| = jS"! = iV hold. 
Next, one calculates the potential as 

V{r)=Eq{kr), (6) 

where 

2d f Kix,x) \ 



X dx 
with 

Kix,y) = Y,AL{x)uUy), (8) 

LGT 

where the coefficient functions Al{x) are calculated by solving the system of linear equations 

E, , , W\ul{x),vi{x^ . X , „ 

with 

/ TTX / TTX 

U-aix)^ J—J^j^^_(x), Vn{x) = J —Y„+l{x), UESUT (10) 

being the regular and irregular solutions of the dimensionless free radial Schrodinger equation for the nth 
partial wave (Riccati-Bessel functions fW). Note that the set T of shifted angular momenta L contains 
real or (in case of complex phase shifts) even complex numbers. 

Using the behaviour of these functions at the origin and infinity, one easily derives the following 
properties of the Cox-Thompson potential. At the origin we get [13], in general, a constant 

q(0) = Q-2il-Q)J2J2^- (11) 

with the matrix Gli = {L — and the constant quantity Q — X^LerX^jes T+K- 
At the infinity the CT potential in general falls off like an inverse power of two 

1 2 

q{x oo) = ^(q!t cos(2a;) — /3t sin(2a;)) -) — ^(aT sin(2a;) -I- /3t cos(2a;) -I- jt) (12) 

where the coefficients depending on the set T are defined by the following relations 

ckT = ^ ^ (flL cos - 6l sin , (13) 

LeT 



LeT 

7T = ^ ^ (^Lcosi^ - aLsini^) , (15) 



LeT 
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with the coefficients aL,bL defined via the asymptotic expansion functions A'^lx) — Al(x — > oo) which 
admits a periodic expansion as 



A%{x) = ql cos(2;) + 6i sin(a;) 
because of the form of the asymptotic version of equation (jH]) 



l{l + I) - L{L + I) 



I £ S. 



(16) 



(17) 



3 Extensions of CT method to long ranged potential 

In nuclear physics one frequently encounters charged particle scattering. The appearance of the long 
ranged Coulomb potential requires an appropriate extension of the CT method which will be carried out 
in various ways in this section. All the different extensions rely upon the two-potential formalism |17j. 

Consider a spherical potential V{r) which can be written as a sum of a short-ranged interior (or 
nuclear) part V and a long-ranged exterior (asymptotic or reference) part V^^'^lr), i.e., 



T/(r) =F(r)+F(")(r). 
Then the phase shift itself is also split into two parts 



(0) 



(18) 



(19) 



where ^J"'' is the phase shift due to y^^-' and Si corresponds to the phase shifts caused by V (in the 
presence of V^°^). If, for instance, V{r) = 0, Si — VL 

If the short ranged part of the potential is zero beyond a finite distance Ta, 



Vir) ^0, r>r^ 

then, accordingly, the radial wave function can be written in this region as 



■^i{k;r> ra) - Q ( (fcr) + tan <5zG|"^ (fcr) 



(20) 



(21) 



In the latter equation the functions F^^^^ and cf^ mean, respectively, the regular and irregular solutions 
of the dimensionless radial Schrodinger equation ([2]) with potential q{x) — E'^^V''^\xk^^). 

3.1 Phase transformation method (PCT) 

We may design a modified potential V{r) which is identical with the original one (up to a constant energy 
shift) within the interior region < r < ra and zero outside 



Vir) 



V{r)-V{r,), r<r,, 
0, r > ra. 



(22) 



Thus if equation (pn| applies the V interior potential can be deduced from V. The radial scattering wave 
function of this problem at the shifted energy Eb — E ~ V{rii) = h^k^/2m is identical with the original 
one in the inside region 

^i(fcB;r-) = *i(fc;0, r<r^ 



but differs in the outside region 

^i{kB;r) = Ci (ui{kBr) - tan(5;u;(fcBr)^ 
The equality of the logarithmic derivatives at r = ra 



r > r„ 



1 ~ 



— log ( -*;(fcB;r 



1 



-log [-^i{k;r 



(23) 
(24) 

(25) 
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with ^';(fc;r) from ((2T|) gives the transformed phase shifts 5i{kB) which can be used to (re)construct 
the short ranged potential V{r) ([21]). To get the whole original potential one simply adds the reference 
potential 

if it is known. 

The phase shift transformation method enables the use of the CT method as formulated in the 
preceding section. It requires that the modified potential V{r) be constructed from the set {Si} by the 
actual inverse method as accurately as possible. But no exact fixed energy inverse scattering method has 
yet been formulated which is able to produce a short range potential that is exactly zero beyond a finite 
radius, if a finite set of phase shifts is employed. Therefore, one always introduces a small error into the 
potential (re) construction when applying the phase shift transformation procedure. 

The phase shift transformation method has been introduced by May, Miinchow and Scheid [B] and 
is widely applied in case of the modified Newton Sabatier method. It is applied also to the CT inverse 
procedure to treat Coulomb scattering and we shall call this combined method as the PCT method. 



3.2 Generalized CT scheme 

Based on the two potential formalism one can derive a generalized CT (gCT) scheme which employs the 
given nuclear phase shifts {5i} for constructing the short-ranged (nuclear) potential V{r) but avoids use 
of a matching radius ra. The derivation starts with the ansatz for the input symmetrical kernel of the 
Gel'fand-Levitan-type integral equation 



g(x,x')=^7^^r(x<)G|")(a;>), (27) 
les 



where a;<(a:>) denotes the lesser (greater) of x,a;'. Then, by proceeding through the usual steps P^fTl] . 
one arrives at a system of nonlinear equations identical in structure to ([3]) 

^"'' = T^ or tan(JO= ,^^f+% I ^ S, (28) 

1 — iK-i 2 + i(K.i — K-i ) 



with 



and 



LeTA'es 



1 1 I sin (?-i)| + 5i°)-5(°) 



^COS 



L{L + 1) - l{l + 1) I cos [{1 - L)f + - Sl'^) 



I eS, LeT (30) 



for the set T of the shifted angular momenta L where the relations S OT — and |r| = \S\ hold. From 
the set T one calculates the nuclear potential as 



, , 2 d f K(x, x) \ 



V{r) = Eq{kr), (31) 

where 

X dx 
with 

if(:r,y) = ^Ai(x)Ff (y). (33) 

The coefficient function Al {x) is calculated by solving the system of linear equations 

W[Fj^\x),G['\x)] 
1(1 + 1) - L{L + 1) 



From this gCT scheme several developments are possible as the actual form of the reference potential 
has not yet been specified. 
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3.3 Coulomb reference potential method (CCT) 

If one sets the reference potential to be the bare Coulomb potential 

V^^\r)=E'^^-^'-^^ (35) 
kr 47r£o r 

with rj being the Sommerfeld parameter then one arrives at the Coulomb CT (CCT) method. In this 
case the regular and irregular reference functions are the regular and irregular Coulomb functions, the 
reference phase shift becomes the Coulomb phase 



(7i ~ — In 
2i 



r(/ + 1 + Z7?) 
r(/ + 1 - 17]) 



I e S. (36) 



Here it should be mentioned that in the course of application of CCT method it may become necessary 
to know the regular and irregular Coulomb functions and Coulomb phases for complex orders. For 
example, in case of non-elastic scattering the phases are complex and therefore, as noted before, the 
L numbers are also complex valued. The Coulomb functions are well-defined for complex orders and 
similarly to the real order case they can be given as power series (for details see e.g. [111 HO] and 
Appendix) . 

It is interesting that contrary to the fact that the asymptotic form of the reference functions contains 
the well known logarithmic term ~riln2kr in the argument, it does not appear in the gCT formulas 
because of cancelation. 

Using the known power series of the Coulomb functions it can be shown that the CCT method gives a 
potential which is proportional to the Coulomb potential near the origin. For one term, i.e. |r| = \S\ = 1, 
we get 



V{r fiiO)^E 

For large r we get 



L{1 + 1) 
1(1 + 1) 



g + 0(l). (37) 



V{r^^) = E^-^Y. E[^^cos]l; cos(ei(fcr) + e,(fcr)) + O I , (38) 
[kr) i^fzT i^s \{'^''^) / 

with 0l(cc) = X — r]lin.2x — L-| + (Jl- 

We see that the CCT method generates an inverse potential that gives a Coulomb-like singularity at 
the origin and produces a damped oscillation around the Coulomb tail at large distances. It is free of the 
matching parameter r^ and requires just the nuclear phase shifts 5i which are derived by the usual phase 
shift analysis procedures. 

3.4 Modified Coulomb reference potential method (MCT) 

In order to obtain an inverse potential that is finite at the origin, instead of being singular there, we can 
modify the Coulomb reference potential accordingly. This reference potential is constant in the interior 
domain and purely Coulombic outside. This modified Coulomb potential is the same as that employed 
by the phase transformation method for Coulomb asymptotics and reads as 

V^'\r)^[lp If I- . (39) 



To this reference potential there belong the following regular and irregular radial wave function 

AO) 



Fr{kr)={ "'(\A^-fc'')' -<'^- 

ap(o)Fi{kr) + l3p(o)Gi{kr), r > r^,, 



G\°\kr)^{ "'VV^ /cr. --y ■ (4^) 

a^mFiikr) + l3^(o)Gi{kr), r > r^,, 



and reference phase shift 



5l°^ =ai + arctan ( \ . (42) 
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The coefficients a^io) , Ppio) , a^co) , P^m can be calculated from the equality of the inner and outer wave 
functions and their derivatives at the matching radius ra. For example, the two coefficients necessary for 
calculating the reference phase shift are 



-^Fr F(x^) _ Fljx^) ' ^^_F,^"' Gi(x^) _ G;(xJ 

Gi{x^) G[{x^) Fi(x^) f;{x^) 

with = fc?'a and prime denotes derivation with respect to argument. 
The total phase shifts can be written in two different ways 

d,=S,+a,=Sf'^^ + Sl°^ (44) 

where 5i means the nuclear phase shifts given as data and Jj^'^'^ is to be used to perform the CT inverse 
calculation outlined above. 

The potential obtained by the MCT method has a finite value at the origin. Because this method 
employs a similar reference potential as the PCT method, the results provided by the two methods should 
also be very similar, although quite different functions are used in the calculations. The advantage of 
the MCT over the PCT lies in that it does not involve the small error in the phase shift reproduction 
inherent to the PCT. 



4 Applications to (re)construct effective potentials 
4.1 Synthetic phase shifts 

To illustrate the general applicability of the long-range CT inversion procedures we shall first use them 
to reconstruct model potentials. We model the a - a scattering with two slightly different potentials: 
the first is finite at the origin and the second is singular (describing a possible non- locality) . 
The first model potential is given by 

Vir) ^Uir) + Vcir), (45) 
where the nuclear interaction is described by a Woods-Saxon form 

U{r) = Uo[l + d'^y\ (46) 
and the Coulomb interaction is represented by the potential of a homogeneous charged sphere as 

Vcir)^\ (47) 
I r>Rc. 

For the various parameters we choose the following values: Zi Z^ — 2, A\ ^ 4, A2 = 4, Uq = —20 
MeV, R = Rc = 2.0636 fm, a = 0.25 fm. 

We have calculated the nuclear phase shifts 61 at energies i? = 25 and 35 MeV (see tables [T] and 
[2]). These phase shifts are then used as input data for the various CT calculations. The results for the 
L- values obtained using the nonlinear solver of MATLAB are shown in tables [T] and [H The corresponding 
potentials are displayed in figure [TJ 

One can observe that as expected the PCT and MCT procedures give almost the same results and the 
CCT inverse potential is divergent at the origin. By increasing the scattering energy more phase shifts 
become available for the inversion (see tables [T] and [5]) and the potential reproduction becomes better. 

The second model is obtained by adding the singular potential term 

^-g(r) - ^ (48) 

to the previous model, i.e., 

V{r) = U{r) + Vc{r) + Fsing(r). (49) 

The reconstruction of phases is listed in table [3] and the potentials are shown in figure [2] We see that the 
CCT potential follows nicely the model potential in the singular domain near the origin while the PCT 
and MCT methods are unable to reproduce the singularity although their phase shift reconstruction is 
good. 
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Table 1: Model data {S°"^ phase shifts), inversion results {L shifted angular momenta and Aj; differences 
between the model phase shifts and the ones given by the various methods|) of the a - a scattering at 
Ec.m. = 25 MeV. The matching parameter used in the PCT and MCT procedures was set to = 10 fm. 



I 




T iVlL/ i 















-1.5622 


-0.9225 


-0.9219 


1.2989 


0.0085 


0.0053 


0.0056 


1 


0.5286 


0.5841 


0.5843 


1.1445 


0.0011 


0.0009 


0.0001 


2 


1.6462 


1.7503 


1.7502 


0.8307 


0.0116 


0.0063 


0.0062 


3 


2.9468 


3.0281 


3.0281 


0.2300 


0.0131 


0.0025 


0.0024 


4 


4.0456 


4.1114 


4.1115 


0.0399 


0.0176 


0.0053 


0.0054 


5 


5.0272 


5.0983 


5.0982 


0.0062 


0.0080 


0.0019 


0.0018 


6 


6.0308 


6.0874 


6.0875 


0.0009 


0.0073 


0.0004 


0.0005 


7 


7.0169 


7.0684 


7.0684 


0.0001 


0.0047 


0.0039 


0.0039 


8 


8.0177 


8.0557 


8.0557 


0.0000 


0.0001 


0.0066 


0.0066 


9 


9.0108 


9.0428 


9.0427 


0.0000 


0.0001 


0.0108 


0.0108 



I Note that for the sake of comparison the phases given by the inverse potentials were calculated by 
cutting-off the non-physical oscillations beyond the matching radius ra used in the PCT procedure. 
Without the cut-off the MCT and CCT potentials reproduce the phase shifts within an error of the 
numerical precision. 

Table 2: Model data and inversion results of the a - a scattering at Ec.m. — 35 MeV. The matching 
parameter used in the PCT and MCT procedures was set to ra = 10 fm. 



I 


^CCT 




^PCT 








A(5PCT 





-1.4866 


-0.8389 


-0.8361 


1.1880 


0.0032 


0.0002 


0.0033 


1 


0.4954 


0.6158 


0.6161 


0.9864 


0.0023 


0.0011 


0.0001 


2 


1.6375 


1.7077 


1.7076 


0.8542 


0.0032 


0.0027 


0.0027 


3 


2.8420 


2.9144 


2.9144 


0.4246 


0.0069 


0.0031 


0.0033 


4 


4.0002 


4.0677 


4.0677 


0.1111 


0.0030 


0.0017 


0.0016 


5 


5.0256 


5.0873 


5.0873 


0.0231 


0.0079 


0.0043 


0.0044 


6 


6.0277 


6.0873 


6.0873 


0.0043 


0.0001 


0.0002 


0.0002 


7 


7.0199 


7.0713 


7.0713 


0.0008 


0.0032 


0.0019 


0.0019 


8 


8.0176 


8.0627 


8.0627 


0.0001 


0.0047 


0.0029 


0.0030 


9 


9.0124 


9.0492 


9.0492 


0.0000 


0.0011 


0.0009 


0.0009 


10 


10.0115 


10.0421 


10.0420 


0.0000 


0.0034 


0.0030 


0.0031 


11 


11.0083 


11.0315 


11.0316 


0.0000 


0.0051 


0.0044 


0.0044 


12 


12.0082 


12.0267 


12.0267 


0.0000 


0.0007 


0.0067 


0.0067 



4.2 Experimental phase shifts 

Our goal is to assess the effective central potential governing the p — a scattering events. Comprehensive 
data of phase shift analysis of the p — a scattering has been presented by Ali, Ahmad and Ferdous in 
pi] . Because of the spin-orbit coupling both spin- up 6'^ and spin-down phase shifts contribute to the 
scattering amplitude at each partial wave. In case of weak spin-orbit coupling the combined phase shifts 

+ + (50) 

are characteristic of the underlying central potential [22!, and we shall use them as input for the CT 
procedures. 

We have inverted the phase shift data of Ali et al. Since the inverse potentials exhibit similar charac- 
teristics we present only some of the inversion results here. 

In figure[3]the inverse potentials yielded by three CT methods at Slab — 17.45 MeV proton energy are 
depicted. The results at this energy are representative of the potentials recovered below the i?iab = 22.94 
MeV, a + p ^ d-l-'^IIe inelastic threshold. As we see the MCT and PCT potentials are almost identical 
and strongly resemble a Woods-Saxon form. The range and strength of all the three potentials are similar 
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r (fm) r (fm) 



(a) Ec.m. = 25 MeV (b) _Bc.m. = 35 MeV 

Figure 1: Model with potential (|45l) and inverse potentials yielded by the CCT, MCT and PCT methods 
(labeled accordingly) at i?c.m. = 25 and 35 MeV center of mass energies. The model potential is non- 
singular at the origin. 



but the CCT potential is different in shape: it possesses a repulsive core. 

Note that a repulsion at small distances can be theoretically accounted for as the manifestation of the 
Coulombic non- locality [2 3) . It can also be seen that this repulsion core stabilizes the inverse potential in 
the sense that the amplitude of the asymptotic oscillations is diminished compared to the non-repulsive 
MCT/PCT results. This is also the reason why the phase shift reproduction with a given precision of 
the CCT potential is better than that of the MCT/PCT potentials (see table [??|) . 

In figure S] numerous other CCT and PCT potentials are shown at various energy values including 
those above inelastic threshold. (The MCT potentials are not shown because they coincide within the 
width of line with the PCT results.) Apart from the Coulombic singularity at the origin, the potentials 
have a similar range of 3 - 4 fm and strength of 50 - 70 MeV (PCT) and 50 - 160 MeV (CCT). The 
imaginary part is much less compared to the real part. The reproduction of phase shifts (not shown) 
gets better at higher energy in case of the PCT potentials because of the fixed cut-off radius which in 
principle does not apply to CCT (and MCT) method. Without use of this radius the CCT (and MCT) 
potentials give back the input phase shifts exactly. 



5 Conclusions, summary 

We have developed the Cox- Thompson (CT) fixed energy inverse quantum scattering method into various 
directions in order to make it appropriate for treating long range interactions. By explicit calculation we 
have shown that the modified (MCT) and phase transformed (PCT) methods yield practically the same 
potentials although quite other functions are involved in the calculation. These potentials are generally 
finite at the origin therefore they can be used for cases with no singularity at small distances. The 
Coulomb (CCT) method produces a characteristic Coulomb singularity at the origin. This method can 
be applied with success when nonlocality or a repulsion hard core plays a role in the interaction. 

We have applied the methods first to model cases with and without singularity. The model potentials 
describing a — a scattering are well reproduced by the new methods according to their characteristic 
properties concerning the finiteness or infiniteness at the origin. Then, the experimental p — a phases of 
Ali et al [21' have been inverted with the result that the CCT potentials showing up a repulsive hard core 
reproduce better the input phase shifts than the PCT potentials which exhibit a Woods-Saxon shape. 
Whether or not the repulsive core is present between the proton and alpha particle when being scattered 
by each other should be clarified also using independent source of information (see e.g. pTll^l^ 1^127) ). 
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Table 3: Model data and inversion results for the divergent model potential at Ec.m. = 25 MeV. The 
matching parameter used in the PCT and MCT procedures was set to ra = 10 fm. 



7 

t 


rCCT 


rMCT 
ij 


rPCT 




A ACCT 


AO; 


A i;PCT 

Ad; 


u 


— 1.5257 


—0.8525 


A O C OO 

—0.8522 


1, 


,2121 


A A A A f 


A AAO O 


A A A A 1 


1 


O.oiib 


U.OODZ 


0.5862 


1. 


,1314 


A AA1 1 
U.UUll 


A AAAO 


A AAAi? 


2 


I.ddOd 


1.7506 


1 A C 

1.7505 


0. 


,8245 


A A A c: 1 

0.0051 


A AA /I O 

0.0043 


A AA A O 

0.0042 


6 


2.9505 


3.0290 


3.0290 


0. 


,2284 


A A A/^ 1 

0.0061 


A AAOO 

0.0028 


A A AOT 

0.0027 


A 

4 


4.U4U1 


4.11(J8 


/111 AA 

4.11(J9 


0. 


,0395 


A AAA A 

0.0094 


A AAf? A 

0.0060 


A AAf? A 

0.0060 


5 


5.0293 


5.0978 


C AATT 

5.0977 


0. 


,0060 


A AAO/^ 

0.0036 


A A A 1 O 

0.0013 


A A A 1 O 

0.0013 


R 
U 


u.uzoo 


U.UO ( u 


R 0871 
U.UO ( 1 


0, 


,0008 


n 0040 


no"?! 

U.UUOl 


no"?! 

U.UUOl 


7 


7.0179 


7.0677 


7.0677 


0, 


,0001 


0.0030 


0.0023 


0.0023 


8 


8.0163 


8.0553 


8.0553 


0, 


,0000 


0.0001 


0.0008 


0.0008 


9 


9.0113 


9.0413 


9.0413 


0, 


,0000 


0.0020 


0.0020 


0.0020 


10 


10.0105 


10.0325 


10.0325 


0, 


,0000 


0.0031 


0.0024 


0.0025 


11 


11.0077 


11.0224 


11.0224 


0, 


,0000 


0.0018 


0.0019 


0.0019 


12 


12.0073 


12.0155 


12.0155 


0, 


,0000 


0.0008 


0.0014 


0.0014 


13 


13.0056 


13.0093 


13.0093 


0, 


,0000 


0.0005 


0.0010 


0.0010 


14 


14.0053 


14.0065 


14.0066 


0, 


,0000 


0.0003 


0.0006 


0.0006 


15 


15.0043 


15.0044 


15.0043 


0, 


,0000 


0.0002 


0.0003 


0.0003 



Table 4: Inversion results of the experimental p — a data at i?iab = 17.45 MeV. The matching parameter 
used in the PCT and MCT procedures was set to Ta = 7 fm from which distance also the MCT and CCT 
potentials have been replaced by the pure Coulombic tail. 



I 




jJACT 


^PCT 




Ad;CCT 




Ad;PCT 





-1.7204 


-1.6703 


-1.6705 


1.7240 


0.0139 


0.0219 


0.0198 


1 


0.5252 


0.4537 


0.4534 


1.4839 


0.0249 


0.0589 


0.0597 


2 


2.0328 


2.0210 


2.0211 


0.0760 


0.0206 


0.0107 


0.0108 


3 


3.0475 


3.0881 


3.0882 


0.0229 


0.0250 


0.0762 


0.0765 



A Coulomb functions of complex order 

In certain applications of both the CCT and MCT formulations of the Cox-Thompson inverse scattering 
method for long-ranged potentials it is necessary to evaluate the regular and irregular Coulomb wave 
functions for complex orders. Based on jl9l 120) the evaluation is accomplished by using the power series 
given below. 

The regular Coulomb wave function for L G C complex order is given by 
P , , 2^d-?[r(L + l + 2>7)r(£ + l-z77)]^ ^ r 

Fl{x) = p^^^ ^ 2) ^ Cjimx^ . (51) 

and the constants are defined by the recursion 

CUi - 1> (52) 

Cl+2 = (53) 

- u'k^^ ^■>^ + 2. (54) 

This production of the Fl{x) is a simple analytic continuation of the formulae in |16) . 

However the irregular Coulomb wave function cannot be given by such a simple generalization. We 
utilize therefore the fact that the Coulomb wave functions and the Whittaker functions are related to 
each other linearly for fixed L and rf. By means of the Whittaker functions the regular Coulomb wave 
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Figure 2: Model with potential (gS]) and inverse potentials yielded by the CCT, MCT and PCT methods 
(labeled accordingly) at E'c.m. = 25 MeV center of mass energy. The model potential is singular at the 
origin. 



function is also given. 



Fl{x) =A{L,r])M{ir],L+l,2ix), (55) 

Gl{x) =iA{L,r])M{ir],L+^,2ix)+B{L,T])W{iT],L+^,2ix) (56) 



with 



A{L, r;) = ^-''^"^'^^'^ [r(^+^g.)r(L+i-..^ (57^ 

The Whittaker functions M(a, 6, y) and W{a, b, y) are well-known for complex arguments (see [TB] chapter 
13). The Gamma function was calculated by a Lanczos series approximation. 

The 77 parameter appearing in the above formulae is the Sommerfeld parameter given with the quan- 
tities discussed earlier as 
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(a) Real part of the CCT potential. 
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(c) Real part of the PCT potential. 
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(b) Imaginary part of the CCT potential. 
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(d) Imaginary part of the PCT potential. 



Figure 4: Complex-valued inverse -p — a potentials yielded by the CCT (a,b) and PCT (c,d) methods at 
various -Eiab proton energies below and above the inelastic threshold using the experimental phase shifts 
of i21j. 
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